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ABSTRACT 

Spatial autocorrelation analysis of 1 2 quantitative characters among 50 populations of Hemernai/lis 
hitkaummis was conducted to understand better their geographic variation patterns. I here are two 
types of relationships between the mean values of the characters measured and their spatial 
autocorrelations: (1) significant heterogeneity of means with significant autocorrelation (mostly floral 
characters); and (2) significant heterogeneity of means with no significant autocorrelation (mostly 
vegetative characters). Perianth size shows significant north-south clinal variation and width of inner 
anti outer perianth show a typical monotonic decline from significant positive autocorrelation ar 
interpopulation distance of b" 7 to 107 km to significant negative autocorrelation at distances or 221 


to 557 km. The results might result from the combinations of gene flow, genetic drift, and/or selec¬ 
tive forces operated in populations of H. hak-uuncnsis. 


RESUMF.N 


Se realize) un analisis tie correlation espacial de 12 caracteres cuantitativos en 30 
poblaciones de Hemerocallis hakianwim para comprender mejor sus patrones de variation 
geografica. Hay dos tipos de relaciones entre los valores medios de los caracteres mediclos y 
sus autocorrelaciones espaciales: (I) heterogeneidad signilicativa de las medias con 
autocorrelation significativa (la mayorfa en caracteres florales); y (2) heterogeneidad 

c J * 

sienificativa de las medias sin autocorrelacion significativa (la mayoria en caracteres 
vegetativos). El tamano del perianto muestra una variacion clinal norte-sur significativa y 
la anchura del perianto interno y externo muestra una disminucion monotonica tipica 
desde una autocorrelacion positiva significativa en distancias interpoblacionales de ()7 a 
107 km hasta una autocorrelacion negativa significativa en distancias de 221 a 357 km. 
Los resultados podrfan derivarse de las combinaciones de 11 ujo genetico, deriva genetica y/ 
o fuerzas selectivas que operan en las poblaciones de H. haknunemis. 
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[INTRODUCTION 


In Korea, five species of Heweroa/His are recognized on the basis of mor¬ 
phological analysis (1 2 quantitative and seven qualitative characters), field¬ 
work, and examination of Japanese herbarium specimens (Chung & Kang 
I 994; Kang & Chung I 997a). These are H. haknutiemh Nakai, Id. tbinibergii 
Baker, //. middendorffii 1 r. et Mey., Id. hongdoeusis M. Chung iN S. Kang, and 
Id. taeanensis S. Kane & M. Chune. 


Hetuet'ouillis hcikmimmis is Korean endemic species that commonly grows 
at the margins of pine-oak forests, grasslands, and disturbed habitats such 
as chest-nut orchards in southern, central, and northwestern Korea (Chung 
(St Kang 1994). Among the four Korean / Ici/uroudlis species, populations 
of II. hiikiumemh are the most variable in their inflorescence and floral 



morphology (Kang & Chung I 994). In addition, Chung and Kang (1 994) 
noted that populations of II. bak/z/ntensis with larger (lowers are generally 
located in the central Korean Peninsula, while populations with smaller 
flowers occur in southern Korea. 1 hese patterns of geographic variation in 
morphological characters of the species require further study for a better 
understanding the species. A detailed study of the patterns of geographic 
variation of quantitative morphology in II. bctkinmams is reported here. 

1 he principal method of analysis for this purpose is spatial autocorrelation 
analysis. Spatial autocorrelation statistics are summary measures of the de¬ 
pendence of the value of a particular variable at one location on the value of 
that same variable at other nearby locations. The most commonly used 
measure is Moran s (1 950) autocorrelation coefficient /, in which geographic 
neighbors are compared in terms of their deviation from the mean of i 
observations. The method has been applied in a variety of fields including 
biogeography, political science, ecology, archaeology, and econometrics (Cliff 
& Orel 1 981; Griffith 1988). Recently, this method of chi iracter analysis 
has been employed in the analysis of morphological characters and the spa¬ 
tial dispersal patterns of genetic variations because it can provide a detailed 
picture of geographic variation in characters (Sokal <N Oden I 978a,b; Sokal 
et al. 1986; Jensen 1986; Chung 1995, 1996). Sokal and Oden (1978a) 

described five correlogram patterns of geographic variation of characters, 
(.line is a gradual change from the highest positive / to the lowest negative 
/, across all distance groups. Depression is nothing but a circular cline, 
with the lowest negative / between the shortest and longest distances, not 
the longest distance. Double depression is on the whole quite similar to 
that of the depression, with low / m a certain distance from which another 
low / was separated. Intrusion pattern shows positive /s at a certain dis¬ 
tance, because of the homogeneity within the distance, with a sharply de- 
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creasing / bordering the distance. Crazy quilt pattern shows high /s sur¬ 
rounded by low Is, vice versa, which yields no significant Is at all. 

The purpose of this study is to determine whether or not the 12 quanti¬ 
tative characters recorded for H. hakuunensis by Chung and Kang (1994) 
and Kang and Chung (1 994) show any geographic patterns such as dines, 
depression, double depression, intrusions, and crazy quilt patterns (Sokal 
& Oden 1928a). 


MATERIALS AND METHODS 


Thirty populations representing the geographic range of H. hakuunensis 
were used for this study (Fig. 1). Measurements of 12 quantitative charac¬ 
ters (Table 1; also see Table 2 in Chung & Kang 1994) were taken on each 
of 10 randomly selected individuals directly from their natural habitats 

from 1991 to 1995. Voucher specimens of all collections are deposited at 

the herbarium of Gyeongsang National University (GNUC). 

For spatial autocorrelation analysis, mean values were assigned to each 
population for the 12 characters. Every possible pair of populations was 
considered as a join or a connection and was assigned to one of six distance 
classes based on the geographic distance between them. These six distance 
classes were constructed by equalizing sample sizes among the classes. The 

distance cl asses are 0<76, _7 G< 1 02, 102<13b, 136<169, 169<221, and 

221 <357 km. Moran’s / values were calculated for interpopulational dis¬ 
tance classes by 


/ = N ZI {W V] I, Zj) (I I 




, J >-' 


• • 


■ ) 


i ) 


(Sokal & Oden 1978a). N is the number of populations, W jj is the join on 
weighting matrix, where W;j is set as one if 7th and /th population are in 
the distance class and zero otherwise, Zj - X\ - X, Zj = Xj - X, the variables 
Xj and Xj are the mean scores for /th and /th population, respectively, and 
X is the mean score for all populations. The value of / ranges between + 1 
(complete positive autocorrelation, i.e., paired populations have identical 
values for all characters) and - I (complete negative autocorrelation). Each / 
value was used to test significant deviations from the expected values, E(/) = - 


1/(N-1) (Cliff & Ord 1981). The neighboring populations in the distance 
class considered, with a significant positive value of Moran s / have similar 
scores, whereas those with a significant negative value should have differ¬ 
ent scores. Overall significance of individual correlograms was tested using 
Bonferroni s criteria (Sakai & Oden 1983)- All calculations and statistical 
analyses were performed using the SAAP program (ver. 4.3) written by D. 
Wartenberg. 
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Table I. List of 12 quantitative characters used in the spatial autocorrelation analysis. 


Acronym 

j 

Character derivation 

Unit or category 

PSH 

Plant (scape) height 

m 

LIF 

Length of inflorescence minus flowers 

cm 

NBS 

Number ot bracts/scape 

# 

NFS 

Number of flowers scape 

# 

LLB 

Length ol the lowest bracts 

cm 

WLB 

Width ol the lowest bracts 

cm 

LPO 

Length of the perianth tube enclosing an ovary 

cm 

LIP 

Length of the inner perianth 

cm 

WIP 

Width of the inner perianth 

cm 

LOP 

Length of the outer perianth 

cm 

WOP 

Width ol the outer perianth 

cm 

WWL 

Width ol the widest leaves 

cm 


RESULTS 


Correlograms for each character are presented in Figures 2—3. Moran’s / 
was significant in 26 of 72 cases (36%). The overall correlogram lor width 
of the lowest bracts (WLB), length of the inner perianth (LIP), width of 
inner perianth (WIP), length of outer perianth (LOP), width of the outer 
perianth (WOP), and width of the widest leaves (WWL) was significant (P 
< 0.01). As presented in Figures 2-3, plant (scape) height (PSH), length ol 
inflorescence minus flowers (LIF), and length of the perianth tube enclos¬ 
ing an ovary (LPO) show depression correlograms (Sokal & Oden 197.3b). 


On the other hand, number ol bracts per scape (NBS), number ol Bower per 


WIP. and WOP represent clinal variation with two, WIP 


WOP. 


that each Boral size character shows a similar pattern of geographic varia¬ 
tion. For example, LIP and LOP, WLP and WOP (fig. 3) revealed significant 
positive autocorrelation (localities connected for that distance class had simi¬ 
lar values) lor distance classes 1 to 3, followed by a steady decline until 
significant negative autocorrelation (localities for that distance class have 


dissimilar values) are recorded lor distance class 3. Although LIP and LOP 
show increase in autocorrelation between distance classes 5 and 6, significant 
negative autocorrelations are also observed on these distance classes (Fig. 
3). As seen in Figures 2-3, length of the lowest bract (LLB) has a crazy 
quilt correlogram pattern with no significant autocorrelations, whereas 
width of the lowest bracts (WLB) and WWL display intrusion patterns. 


DISCUSSION 

There were two types of relationships between means and spatial 
autocorrelations observed. The first pattern of significant heterogeneity of 
means with significant autocorrelation was observed for WLB, LIP, WIP, 
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Fiti. 2. Corri-lograms tor PS I I, LIP, NBS, NFS, LLB, and WLB. SivniLic ant autocorrelation 


coefficients (P < 0.05) are indicated bv asterisks. 


LOR WOP, and WWL. lhc second pattern ol significant heterogeniety of 
means with no significant autocorrelation is found for the other six charac¬ 
ters. Sokal and Oden (1 978b) suggested biological explanations for each of 
the two types of patterns. The first could result from migration, founder 
effects, and selective agents. The second type of relationship between means 
and autocorrelations could result from drift, weak migration and selection, 
frequent local extinctions followed by the establishment of new patches by 
founders, or from selection in isolated, patchy environments. Except for 
LLB, all characters show short distance positive autocorrelations, indicat¬ 
ing that migration could be operating over the distances. Certainly in some 
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Fici. 3- Correlograms forLPO, LIP, WIP, LOP, WOP, and WWL. Significant autocorrelation 
coefficients (P < 0.05) are indicated by asterisks. 


of the populations examined there are small patches whose diameter is less 
than interpopulational distances. The width ol the lowest bracts and the 
longest leaves yielded a relatively long distance positive autocorrelation, 
which indicates these characters occur with circular or symmetrical gradi¬ 
ents. It is of interest to note that the width of inner (WIP) and outer (WOP) 
perianth show a typical monotonic decline from significant positive 
autocorrelation at 67 to 102 Km to a significant negative autocorrelogram 
from 221 to 357 Km. Monotonicity is expected when separation by dis¬ 
tance and gene flow between neighbor populations is the main factor. 













































































































Fhe results of this study support the suggestions by Chung and Kang 


(1994) that there is a gradual transition of the floral morphology of //. 
hcikmniensh from larger flowers in the central Korean Peninsula to smaller 

o 

flowers in the southern Korean Peninsula. A similar pattern of geoclinal 
variation of floral morphology is previously reported in H. cl/nuortieri com¬ 
plex in japan (Noguchi 1986), who noted that the flower tube enclosing 
the ovary and the length of the internal as well as external perianths exhib¬ 
ited a geoclinal variation. For example, northern populations possessed much 
shorter, somewhat thicker flower tubes and smaller perianths in compari¬ 
son with those from southern populations. 

In general, vegetative characters have much showed phenoplasticity than 
do reproductive characters in response to environmental and developmen¬ 
tal factors (Stebbins 1 9A()). It might be more reasonable to discuss the 
biolo gical meanings ot the results using reproductive characters. Accord¬ 
ing to the criteria of Sokal and Oden (I 978b), patterns such as those seen 
in several reproductive characters of H. hitkmmemh might result from the 
combinations of gene flow via pollen and seed dispersal, founder effects, 
genetic drift, selective forces in patchy, isolated environments. More re¬ 
cently, we have analyzed the spatial distribution of genotypes using spatial 
autocorrelation of 28 alleles in populations (20 X 2()-m) of II. hakitmiensis. 
A genetic similarity was shared among individuals within 4.5 m distance, 
whereas an overall genetic dissimilarity among individuals beyond the dis¬ 
tance, indicating that individuals in populations of H. httkiunieusis occur on 
a gradient (Chung unpubl. data). Based on allozyme data obtained from 1 9 
populations of //. hcikiimiensh , Kang and Chung (1 997b) further suggested 
that gene movement and genetic drift may be primary factors for shaping 
population genetic structure among populations of H. hakmniemis. 
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